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The variables were defined as: 
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R„ is the radius of the wall for axisymmetric flow (k=l). 
For rectangular coordinates R,, was ignored (k-0). It 
should be noted that Re and time terms were not actual- 
ly non -dimensional. Re had units of (length) 1 , and time 
had units of length. Also, x and y (and Rj remained 
dimensional. So. each term in each of the above bound- 
ary layer equations has units of (length) 1 . 


SOLUTION ALGORITHM 

The momentum, continuity and energy equations were 
differenced as described by Kwon et al.[l] 

For any scalar quantity <}> : 
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The equations were solved at a given time step know- 
ing the solution at the previous time step. Due to the 
parabolic nature of the boundary layer equations, the 
governing equations were solved at each time step by 
marching from station to station in the direction of the 
core flow. Solutions at a given station were obtained by 
solving the boundary layer equations sequentially. First, 
the momentum equation was solved for the velocity 
component in the core flow (streamwise) direction (u). 
Second, continuity was solved for the velocity compo- 
nent in the y direction (v). Then, the energy equation 
was solved for the temperature (T). If the station had 
not converged, the momentum, continuity, and energy 
equations were solved again for the velocity components 
and temperature, using the flow properties from the 
previous iteration. Convergence at a given station was 
obtained when the streamwise velocity components at 
all grid locations at that station converged. After the 
calculation at the station converged, the algorithm 
marched to the next streamwise station to solve for the 
boundary layer properties. This streamwise marching 
continued throughout the entire domain. Once the entire 
solution was obtained for the current time step, the 
boundary layer properties were used to calculate the 
core flow solution at the next time step. 

Grid spacing perpendicular to the wall (y-direction) 
was based on an exponential function. The grid spacing 
was fine near the wall to better resolve the gradients at 
the wall. The grid spacing was course away from the 
wall where fine resolution of the gradients normal to the 
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core stream were not necessary. 

The derivative of the velocity at the wall (used to 
calculate wall shear stress) was determined using a 
second order approximation. Values at j=l (wall), j=2, 
and j*3 were used. 
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The above algorithm worked well for flows with weak 
viscous/inviscid interaction, since the downstream influ- 
ence could be neglected. However, when the flows 
were strongly interacting, such as those with strongly 
adverse pressure gradients or separation, downstream 
conditions had to be considered. The downstream influ- 
ences were treated in two ways ( I ) the pressure gradient 
was differenced as a weighted average of forward and 
backward differences, and (2) a global pressure sweep 
was used. Both of these techniques were based on the 
method of Davis and Barnett[2], 


bolized Navier-Stokes (PNS) equations to convey up- 
stream influences. They had modified the pressure 
gradient terms using the downstream pressure. Because, 
the downstream pressure had been used, a global pres- 
sure sweep had been needed throughout the subsonic 
regions of the solution domain. 

This pressure sweep technique also worked well with 
the uncoupled, unsteady PNS equations described above. 
The fictitious unsteady term that had been used by 
Davis and Barnett was not necessary in our case since 
the flow being modeled was. itself, unsteady. The 
equation for the pressure sweep was then 

P =P -P old +p old ,|»i 
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This pressure technique permitted the modeling of 
moderately separated flows. Downstream effects could 
be quickly propagated upstream. This method was 
made even faster by applying the pressure sweep only 
from a region of weak interaction (downstream of reat- 
tachment) to just upstream of the point of separation. 


Specifically the pressure gradient was differenced as 
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where e is a weighting parameter of the forward and 
backward differencing. This term was required to re- 
move the ellipticity in the PNS equations in strongly 
interacting flows. The quantity e determines what frac- 
tion of the forward difference of the pressure gradient 
that can be included so that the equations remain non- 
elliptic. If the flow was supersonic at a given j location, 
then only backward differencing was used for the pres- 
sure gradient (6=1). If the flow was subsonic then at a 
given j location then the following expression for e was 
used. 


In addition, flaring of the advective term parallel to 
the core flow was used when the flow was separated. 
This was handled by taking the absolute value of the 
coefficient in front of the derivative of the advective 
term (RUB) and multiplying it by .1 to make it smaller. 
Again this operation was only performed when the flow 
was separated. 

A modified Baldwin-Lomax model was used to ac- 
count for turbulence. This model was a zero equation, 
eddy viscosity model. This model was faster than other 
turbulence models, such as the k-6 two equation model. 
A modified model of Visbal and Knight was used to 
handle the effects of separation. Also to account for up- 
stream turbulence history effects, a relaxation method of 
Shang and Hankey was employed. A detailed descrip- 
tion of the turbulence model used is given in Sakowski 
et. al [3], 
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Most boundary layer methods use an inverse tech- 
nique when the flow is separated. This requires a lot of 
computer time because the core flow and boundary layer 
regions must be modeled simultaneously. A method to 
convey the downstream influences upstream, while still 
maintaining the streamwise marching technique was 
considered. This made the algorithm much faster than 
other separated flow algorithms. Davis and Bamett[2] 
had modified the pressure terms in the supersonic para- 


This uncoupled boundary layer scheme worked well 
when used with the core flow solution already known. 
With this method, we can quickly determine the effects 
of boundary layers within inlets of supersonic jet en- 
gines. Thus, we can model several different geometric 
configurations and flow conditions in a reasonable 
amount of time. 


RESULTS AND DISCUSSION 
This uncoupled technique was compared with a tech- 
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nique in which the boundary layer equations were 
solved simultaneously (coupled). The boundary layer 
equations and interactions with the inviscid core flow 
were set up the same way for both techniques. Each 
technique was executed until it converged to the same 
degree of accuracy. The uncoupled technique converged 
quickly, in 2 or 3 iterations, at a each streamwise sta- 
tion. The uncoupled technique tended to be 2.4 to 4.0 
times faster than the coupled technique. 

For the rest of this section we will attempt to verify 
the boundary layer algorithm and to judge the relative 
importance of various aspects of the algorithm. To do 
this we will look at two examples. The first example 
involves a case where the core flow is prescribed and 
the boundary layer is not interactive with the core flow. 
Specifically, the boundary layer algorithm will be used 
to model the flow conditions measured in experiments 
by Lewis, et al,[4], For this example the core flow 
Mach number distribution was prescribed and was not 
interactive with the boundary layer code. The core flow 
Mach number distribution is shown on figure 2b. Cal- 
culations were started using a zero pressure gradient 
profile with a Reynolds number, based on momentum 
thickness, of about 4800 at x=l 1.5 inches downstream 
of the leading edge. 

A second example involved a model of an actual 
supersonic inlet. In this example the boundary layer 
algorithm was used with the Large Perturbation Inlet 
code (LAPIN). LAPIN was used to calculate the core 
flow properties upstream and within a supersonic inlet. 

In this case the boundary layer code and LAPIN were 
interactive. Figure 1 shows the geometry of the NASA 
- LeRC 40-60 inlet. A free stream Mach number of 2.5 
was used. The predictions of the combined boundary 
layer/core flow code were compared with the experi- 
mental results of Cubbison, et al. 

Dependence on grid spacing 

The first example, with the prescribed, noil-interactive 
core flow, was modeled using various numbers of grid 
points perpendicular to the core flow. If the number of 
grid points perpendicular to the core flow direction was 
too small (less than 20 in our case), the algorithm was 
unstable. Increasing the number of grid points made the 
algorithm stable. Good accuracy, however, is not ob- 
tained until an optimum number of grid points are used. 
When the number of grid points is too low the algorithm 
predicts friction coefficients that are lower than experi- 
mental results. To a point the predicted friction coeffi- 
cient increased as the number of grid points was in- 
creased. However, the friction coefficient does not 


increase monotonically with the number of grid points. 
Beyond a certain number of grid points, the friction 
coefficient actually goes down as the number of grid 
points was increased. The point with the highest fric- 
tion coefficient also seems to be closest to the experi- 
mental observations. This is the optimal number of grid 
points. 

Figure 2a shows the predicted friction coefficient ver- 
sus x for various numbers of grid points perpendicular 
to the core flow. In addition, the predicted results were 
compared with experimental results of Lewis, et al. [4] 
and the analytical model of Cebeci and Smith[6]. In 
this case, if the number of giid points was less than 100. 
the friction coefficient increased by increasing the grid 
points. Beyond 100 grid points the friction coefficient 
dropped by increasing the number of grid points. 100 
grid points was the optimum number in terms of the 
accuracy of the results. However, any number of grid 
points between 50 and 100 yielded reasonable results. 
Obviously, increasing the number of grid points from 50 
to 100 resulted in much longer computing times. The 
best number of grid points to use depends on the degree 
of accuracy that is required and the computing time that 
can be allocated. For instance, if one were interested in 
detailed wall effects, heat transfer or skin friction, it 
may be worth while to use 100 grid points. The accura- 
cy is also dependent upon parameters used in the turbu- 
lence model. The results illustrated in figure 2a are for 
a given set of turbulent parameters used in the Baldwin- 
Lomax model. Sakowski, et al.[3] showed that for this 
particular case, increasing the number of grid points 
perpendicular to the wall required some adjustment to 
the turbulent parameters to yield results which more 
closely matched the experimental data. 

However, if large scale flow properties are of concern, 
core flow Mach number and pressure, 50 points or less 
may be used. This was evident from the results of the 
second example. Figure 3 shows the dependency of the 
number of grid points perpendicular to the wall on the 
predicted core flow in the NASA - LeRC 40-60 inlet. 
Note there are two boundary layers in this case, one on 
the centerbody and one on the cowl. The number of 
grid points refers to each boundary layer, not the total 
for both. As in the first example, when the number of 
grid points is too low (20 or less) the algorithm is unsta- 
ble. in this case, the core flow and boundary layer 
codes were interactive, so the boundary layer results 
affected the core flow properties. For instance, when 
the number of grid points was increased to around 30 
the algorithm was stable but the shock wave was pre- 
dicted downstream of the experimental observation. The 
predicted pressures away from the shock were reason- 
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able even when only 30 grid points were used. When 
the number of grid points was increased to 40 the pre- 
dicted shock wave position was further upstream and 
matched the experimental results better. When the 
number of grid points was increased beyond 40 the 
change in the core flow pressure was negligible. Thus, 
in this case there was no need to use more than 40 grid 
points perpendicular to the core flow. 

Dependence on Convergence criterion 

As described above the algorithm iterated until the 
velocities parallel to the core flow (u> converged. The 
convergence criterion at a given x location was a frac- 
tion of the local core flow velocity (u e ). Obviously, 
tightening the convergence criterion increased the com- 
puting time as well as increased the accuracy. It was 
found that if the convergence criterion was too loose, 
the solution was simply chaotic and unrealistic. Alter- 
ing the convergence criterion had the biggest impact on 
the wall values, such as heat transfer, skin friction pre- 
dictions and momentum thickness. As long as the con- 
vergence criterion was strict enough, the predicted flow 
properties were negligibly affected by the grid spacing 
in the x direction. And. if the core flow properties were 
of concern, the convergence criterion itself had only a 
small effect. In fact the convergence criterion can be as 
high as 10% of the local core flow velocity with only 
slight impact on the core flow properties. Figure 4 
shows the predicted pressure as a function of position 
for flow in the 40-60 inlet for various levels of con- 
vergence. As shown on figure 4 the convergence crite- 
rion only affected the predicted core flow pressure near 
the throat of the inlet. 

Effect of upstream propagation 

In prescribing an inviscid core flow solution, the 
backward pressure sweep was not very significant in 
strongly interacting regions, i.e. shocks and adverse 
pressure gradients. The more significant mode of propa- 
gating effects upstream is the calculation of the pressure 
gradient as described in equation 16. The forward dif- 
ferencing was necessary in adjusting the pressure gradi- 
ent to provide a path for upstream propagation when the 
flow is subsonic. However, only a portion of the for- 
ward difference had be used so that the equations did 
not become elliptic. When the weighted for- 
ward/backward differencing of the pressure gradient was 
omitted, the results deviated from the experimental 
results. 

However, it was found that the results were the same 
if pressure sweep was included or omitted. This is most 


likely due to the imposed core flow solution which 
already included downstream effects when the core flow 
was subsonic. 

For the model of the 40-60 inlet, the pressure sweep 
had no influence on the predicted results. However the 
forward differencing of the pressure gradient had an 
effect. Figure 5 shows the predicted pressure distribu- 
tion in the 40-60 inlet using the pressure gradient de- 
scribed in equations 16 and 17 and the predicted pres- 
sure distribution when only backward differencing of the 
pressure gradient was used ( 6=1 ). As expected, there 
was almost no difference in two results when the core 
flow was supersonic, since e=I was used in supersonic 
regions of the boundary layer in both cases. However, 
when the core flow was subsonic the predictions di- 
verged from the experimental observation when only 
backward differencing was used. 

CONCLUSIONS 

Using the uncoupled boundary layer scheme was 2.4 
to 4.0 times faster than a boundary layer scheme that 
solves the governing equations simultaneously. Thus, 
this method was a quick way to model boundary laver 
effects when a core flow algorithm was available. 

Significant considerations were the grid spacing in the 
direction perpendicular to the wall, the differencing of 
the pressure gradient, and the convergence criterion. 
Changing the grid spacing and the convergence criterion 
had the largest impact when specific boundary layer 
properties were of concern, such and friction coefficient 
and momentum thickness. The Baldwin-Lomax turbu- 
lence model, which was used to model turbulence, was 
shown by Sakowski. et al. [3] to be sensitive to the grid 
spacing, which in turn greatly affected the skin friction 
coefficient. If the convergence criteria was not small 
enough, the boundary layer algorithm was not robust in 
predicting boundary layer properties such as skin fric- 
tion and momentum thickness. These parameters be- 
came highly sensitive to the grid spacing in the stream- 
wise direction when the convergence criteria was not 
small enough. 

However the grid spacing and convergence criterion 
had less of an impact on the overall core flow properties 
than the boundary layer properties. In fact, as long as 
the algorithm was stable, the number of grid points and 
convergence criterion could have many different values 
without having a significant effect on the core flow. 

The differencing of the pressure gradient seemed to 
have the largest impact on the core flow solution. The 
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weighted average of forward and backward differencing 
resulted in predictions closer to the experimental obser- 
vation than when only backward differencing was used. 
Finally, the changes made to the Baldwin-Lomax Tur- 
bulence model, described by Sakowski, et al.[3], were 
significant to the accuracy of the results. 
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Figure 1 Geometry of NASA-LeRC 

40-60 supersonic inlet. 
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Figure 2a 
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Figure 3 Predictions of pressure distribution for Figure 5 

NASA-LeRC 40-60 inJet varying grid 
points perpendicular to wall. 
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Figure 4 Predictions of pressure distributions 

of the NASA-LeRC 40-60 inlet 
varying the level of convergence. 
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Predictions of pressure distribution in 
the NASA-LeRC 40-60 iniet varying 
the pressure gradient parameter <e>. 


7 






REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, .nclud.ng the time tor rev.ewtn g instructions searching ex.stino data sources 
ga henng and maintaining the data needed, and completing and revtew.ng the collection of information. Send comments regarding this burden estimate as^fof thlk 

collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Directorate for information OoeraS^ 

Davis Highway, Surte 1204. Arlington, VA 22202-4302. and to the Office of Management and Budget. Paperwork ReSSTp“ Jefferson 

1. AGENCY USE ONLY {Leave blank) T REPORT DATE V 3 . REPORT TYPE AND DATES COVERED 

July 1992 Technical Memorandum 

4. TITLE AND SUBTITLE ' | 5 FUNDING NUMBERS 

A Fast, Uncoupled, Compressible, Two-Dimensional, Unsteady Boundary 
Layer Algorithm With Separation for Engine Inlets 


6. AUTHOR (S) 

Robert L. Roach, Chris Nelson, Barbara Sakowski, 
Douglas Darling, and Allan G. van de Wall 


WU -505-62-20 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-7065 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 105686 
AIAA-92-3082 


11. SUPPLEMENTARY NOTES " 

Prepared for the 28th Joint Propulsion Conference and Exhibit cosponsored by AIAA, SAE, ASME, and ASEE, Nashville, Tennessee, 

July 6-8, 1992. Robert L. Roach and Chris Nelson, Georgia Institute of Technology, Atlanta, Georgia, 30332; Barbara Sakowski and Douglas Darling, 

NASA Lewis Research Center, Cleveland, Ohio; and Allan G. van dc Wall, Case Western Reserve University, Cleveland, Ohio 44106. Responsible person, 
Douglas Darling, (2 1 6) 433-2370. 

^2aT DISTRIBUTION/AVAILABILITY STATEMENT _ 12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Categories 07, 34, and 61 

13. ABSTRACT (Maximum 200 words) 

A finite difference boundary layer algorithm was developed to model viscous effects when an inviscid core flow 
solution is given. This algorithm solved each boundary layer equation separately, then iterated to find a solution. 
Solving the boundary layer equations sequentially, in this manner, was 2.4 to 4.0 times faster than solving the bound- 
ary layer equations simultaneously. This algorithm used a modified Baldwin-Lomax turbulence model, a weighted 
average of forward and backward differencing of the pressure gradient, and a backward sweep of the pressure. With 
these modifications, the boundary layer algorithm was able to model flows with and without separation. The number 
of grid points used in the boundary layer algorithm affected the stability of the algorithm as well as the accuracy of 
the predictions of friction coefficients and momentum thicknesses. Results of this boundary layer algorithm com- 
pared well with experimental observations of friction coefficients and momentum thicknesses. In addition, when used 
interactively with an inviscid flow algorithm, this boundary layer algorithm corrected for viscous effects to give a 
good match with experimental observations for pressures in a supersonic inlet. 


14. SUBJECT TERMS 

Boundary layer flow; Numerical analysis compressible flow; Intake systems 


15. NUMBER OF PAGES 

10 

16. PRICE CODE 

A02 


17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19. SECURITY CLASSIFICATION 

20. LIMITATION OF ABSTRACT 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 

Unclassified 

Unclassified 

Unclassified 



NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 









